On X-ray Optical Depth in the Coronae of Active Stars 

Paola Testa^, Jeremy J. Drake^, Giovanni Peres^, David P. Huenemoerder^ 

ABSTRACT 

We have investigated the optical thickness of the coronal plasma through 
the analysis of high-resolution X-ray spectra of a large sample of active stars 
observed with the High Energy Transmission Grating Spectrometer on Chandra. 
In particular, we probed for the presence of significant resonant scattering in the 
strong Lyman series lines arising from hydrogen-like oxygen and neon ions. The 
active RS CVn-type binaries II Peg and IM Peg and the single M dwarf EV Lac 
show significant optical depth. For these active coronae, the Lja/Lj/3 ratios 
are significantly depleted as compared with theoretical predictions and with the 
same ratios observed in similar active stars. Interpreting these decrements in 
terms of resonance scattering of line photons out of the line-of-sight, we are able 
to derive an estimate for the typical size of coronal structures, and from these 
we also derive estimates of coronal filling factors. For all three sources we find 
that the both the photon path length as a fraction of the stellar radius, and the 
implied surface filling factors are very small and amount to a few percent at most. 
The measured Lja/LjjS ratios are in good agreement with APED theoretical 
predictions, thus indicating negligible optical depth, for the other sources in our 
sample. We discuss the implications for coronal structuring and heating flux 
requirements. For the stellar sample as a whole, the data suggest increasing 
quenching of Lya relative to Ly/? as function of both Lx/i^boi and the density- 
sensitive Mgxi forbidden to intercombination line ratio, as might generally be 
expected. 
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Introduction 



A fundamental issue in the physics of stellar outer atmospheres concerns the relationship 
between magnetic activity on stars with a wide range of physical parameters and solar 
magnetic activity. How directly and how far does the solar analogy apply to other stars, 
and do any of the underlying physica l processes differ? The X-ray luminosities of late-type 
stars can span several decades (e.g., IVaiana et al.lll98ll ). and these hot coronae are found 
on such a wide range of spectral types that the extrapolation of the now well-studied solar 
corona to the extremes of stellar activity is by no means obvious and could be inappropriate. 
Indeed, "scaled up Sun" scenarios, in which a stellar surface is covered with bright solar-like 
active regions, on ly realise X-ray luminosities 100 times that of the typical active Sun (e.g., 
Drake et al.ll2000l ). The most active stars, with X-ray luminosities of up to 10,000 times the 
solar X-ray luminosity, must have coronae which are structured differently in some way. 

Since coronal structures can be imaged presently only on the Sun, the structuring of 
other stellar coronae is generally investigated through the application of techniques such 



as the study of lightcurves during flares (e.g., ISchmi tt fc Favata 1999 : Favata et aL 



Maggio et al.ll2000 



2001 



Marino et al. 



2003 



2000a 



Reale et al.ll2004j: iTesta et al...20 07'). rotational modulation (e.g.. lBrickhouse et al 



Huenemoerder et al.l 120061). and study of de nsity properties to 



gether with info rmation on emission measure (e.g., iTesta et al.l l2004al . hereafter TDP04; 



Ness et al.ll2004j ). Most of these analyses indicate that the emittin g plasma is rather com- 



pact (scale height < . 5-R^) and localized at high latitude (see e.g.. ISchmitt fc Favatalll999 



Brickhouse et allboOll : iTesta et al.ll2004bl . TDP04); however, the presence of extended coro- 
nal plasma has also been claimed on some stars based on UV and X-ray Doppler studies 
(e.g. 



Chung et al.ll2004J : iRedfleld et al.l 120031 ) 



The search for signs of quenching in strong lines through resonance scattering represents 
a further technique that offers a potentially powerful diagnostic of the sizes of X-ray emitting 
regions; the escape probability of a photon emitted by a resonance line in a low density 
homogeneous plasma is in fact dependent on the line-of-sight path length through the plasma 
region. Signiflcant scattering optical depth can be combined with density measurements to 
obtain an estimate of photon path length within the emitting plasma. 

Several existing studies have explored optical depths of both solar and stellar coronal 
emission lines. Studies of solar X-ray spectra have aimed at probing the optical depth in 
the strong {gf = 2.66) 2p^3d^Pi — 2p^ ^Sq resonance line of Fexvil at 15.01 A as compared 
to nearby weaker Fe xvil lines, t hough with controversial r e sults concerning whether optical 



depth effects were se en or not (IPhillips et al.l Il996l . 119971 : ISchmelz et al.l 119971 : ISaba et al. 
19991 ). In particular, ISaba et al.l ( 19991 ) review recent observational flndings on the opacity 



inferred from the study of the bright iron resonance line at 15.01 A and on the center- 
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to-limb behaviour. Among other issues, ISaba et aLl (119991 ) address the discrepancy they 
find i n the derived direct ion and magnitude of the center-t o-hmb trend (a l so in agreement 
with ISchmelz et aLl 119971 ) , as compared to the findings of iPhilhps et aLl (jl996l ) who find 
that the effect of resonant scattering is decreasing from the disk center toward the solar 
limb, a trend i rrecon ci lable and totally oppos i te to that found by ISaba et al.! (119991 ) and 
Schmelz et aLl (119971 ). iBrickhouse fc Schmela (120061 ) have recently reanalized solar X-ray 
spectra and suggest that previously ignored blends might explain the departure of measured 
ratios from theoretical calculations. 



Resonant scatter ing in stellar coronae has been investigated by jPhillips et al.l (120011 ) 
and iNess et al.l (120031 ) through the analysis of the same transition observed at high resolu- 
tion (A/AA up to ~ 1000) by Chandra and XMM-Newton. Both stellar studies of Fexvil 
transitions fail to find evidence for significant deviation from the optically thin regime, and 
in particular the large survey of stellar spectra analyzed by iNess et al.l (120031 ) show that 
no firm results can be obtained from Fe lines. One exception is the suggestion of resonance 
scattering in Fexvil 15.01 A on the basis of t he observed var iability of line ratios seen during 
a fiare on AB Dor by iMatranga et al.l (120051) . We note that IPhillips et al.l (1200 ll ) attempted 
to derive constraints on emitting region size based on their upper limit to optical depth in 
the coronae of Capella, though, as we discuss in this paper (^, such upper limits cannot 
be reliably used in this fashion because scattering into the line of sight renders optical depth 
measurements themselves only lower limits to the true scattering depth. 



In a previous Letter (iTesta et al.ll2004bl . hereafter Paper 1) we presented results obtained 
using a different approach to the study of coronal optical depth, through the analysis of 
Ne and O Lya to Ly/5 line strength ratios as observed by the High Energy Transmission 
Grating (HETG) on board the Chandra X-ray Observatory. Significant depletion of Lya 
lines to resonance scattering were seen in the spectra of the RS CVn-type binaries II Peg 
and IM Peg. In this paper we follow up on that exploratory study and examine the Ne 
and O Lyman series lines in a large number of active stars (namely, the same sample for 
which the plasma density was analysed in TDP04) in order to survey the X-ray optical depth 
properties of active stellar coronae. 

We discuss in ^the advantages of the Lyman series analysis with respect to the "stan- 
dard" approach using Fexvil lines. The observations are briefiy described in ^ Our 
techniques of line fiux measurement and spectral analysis are described in ^ The results 
are presented in ^ We combine the results of this study with our earlier density estimates 
and discuss these in the context of coronal structure on active stars in §S1 we draw our 
conclusions in ^ 
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2. Resonant scattering in Ne and O Lyman series lines. 



The Fe xvil soft X-ray complex at ~ IsA has been a primary tool to probe coronal op- 
tical depth because of the large oscillator strength of the ^ Sv, -2jr'?)d} 15.01 A reson ance 
line and its consequent prominence in solar spectra. However, iDoron fc Beharl (120021 ) and 
Gul ( l2003bl ) have recently shown that the indirect processes of radiative recombination, di- 
electronic recombination, and resonance excitation involving the neighbouring charge states 
are important for understanding the relative strengths of Fexvil-xx lines. There is thus 
still some considerable difficulty i n reco nciling theoretical and observed line strength ratios. 
Recently, iBrickhouse fc Schmelzl (120061) have found good a greement of solar observed ratios 
with new theore tical calculations (Chen fc Pradhan 2005h. and suggest t hat center-to-limb 



Schmelz et al. 



199 



3 



Saba et al. 



19991 ) are due to chance 



observed trends ( Phillips et al. 19961 : 
rather than to optical depth effects. 

An additional problem in using Fexvil lines as diagnostics of optical depth is that 



Huenemoerder et al. 


2001; 


Drake 


2003: 


Audard et al. 


2003) 



as compared with a solar or local cosmic composition. The ratio of the line-center optical 
depths, Tj/xj, of two lines i and j is given by 



fi\i(j)iAiJm. 



where / is the oscillator strength, A is the wavelength, (p is the fractional population of the 
ion in question, A is the element abundance and m the ion mass. The line optical depth is 
directly proportional to abundance and, in the case of stellar coronae, where only the very 
strongest spectral lines might be expected to undergo any significant resonance scattering, 
any abundance depletions also reduce the sensitivity of lines as optical depth in dicators. For 



coronal abundance s typically found in RS CVn systems (see e.g. reviews by iDrakd 12003 



Audard et al.ll2003l ). we expect resonant lines from the more abundant ions like oxygen and 
neon to be more sensitive to resonant scattering processes than the resonant Fexvil line. 
This is illustrated in Figure [H where we show the relative optical depths for the O vill and 
Nex Lya lines and th e Fexvil ~ 15.01 A res onance line for both a representative solar 



chemical composition (iGrevesse fc Sauvall Il998l ) and for a chemical composition typically 
found for active stars. For the latter, we assumed Ne and Fe abundances 0.3 dex higher and 
0.5 dex lower, respectively, than the corresponding solar values; we note that in several active 
coronae abundance anomalies ev en more pronounced have been found (e.g 
200ll : iHuenemoerder et aDboOlh . 



Brinkman et al. 
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Opticol depth comporison for Fe XVII, Ne X, VIII 




5.0 5.5 6.0 6.5 7.0 7.5 
Log T [K] 



Fig. 1. — The relative sensitivity of Lya lines of Nex, O vill, and the strong Fexvil 15.01 A 
transition to resonant scattering, as a function of plasma temperature. Relative line optical 
depths, normalised to the maximum depth in O Vlll, are shown for two sets of abundances- 
Dashe d lines correspond to the solar photospheric abundance mixture of iGrevesse fc Sauval 
( 1l998l ): filled profiles correspond to coronal abundances typically found in RS CVn's (see 
text). Because of the typically higher Ne abundance and lower Fe abundance found in the 
coronae of RS CVn systems, the expected resonant scattering effects are greater for O vill 
and Nex Lya than for the Fexvil. 
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As discussed in Paper I, the effect of resonance scattering of Lya and Ly/3 pliotons can 
be diagnosed by comparison of tlie measured Lya/Ly/S ratio witfi respect to the theoretical 
ratio. In principle, both n = 2 1 Lja and n = 3 —>■ 1 Ly/5 lines can be affected by resonant 
scattering: when a large optical depth is reached in Lya, an enhanced population of the n = 2 
level can lead to a potentially confusing enhancement in Ly/3 (and Baa) through collisional 
excitation of the n = 2 ^ 3 transition. However, in the limit of fairly small scattering optical 
depth that we expect to characterize stellar coronae, Ly/3 essentially remains optically thin 
and should make a reliable comparison with which to diagnose optical depth in Lyct. 



3. Observations 



Chandra High Energy Transmission Grating Spectrometer (see ICanizares et al.l |2000| . 



20051 for a description of the instrumentation) observations of 22 cool stars covering a wide 
range activity level were analysed. The stellar parameters and particulars of the observa- 
tions were discussed in detail in a companion paper that addressed the coronal densities of 
active stars (TDP04). To the sample analysed in TDP04 we added a set of observations of 
IM Peg obtained subsequent to the first three segments analysed in Paper I and TDP04; the 
parameters of these additional observations are listed in Table [TJ 
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Table 1. Parameters of the HETG observations of IM Peg. 



Obs ID 


Start date and time 


^exp 










[ks] 


[erg s"^] 


[erg s"^] 


2527 


2002-07-01, 15:39:08 


24.6 


2.95x10^1 


3.04x10^1 


2528 


2002-07-08, 02:07:29 


24.8 


2.37x10^1 


2.55x10^1 


2529 


2002-07-13, 06:27:34 


24.8 


2.03x10^1 


2.19x1031 


2530 


2002-07-18, 21:13:58 


23.9 


1.96x10^1 


2.15x1031 


2531 


2002-07-25, 11:30:11 


23.9 


1.99x10^1 


2.20x10^1 


2532 


2002-07-31, 20:45:59 


22.5 


3.00x10^1 


3.23x10^1 


2533 


2002-08-08, 00:50:22 


23.7 


2.57x10^1 


2.79x10^1 


2534 


2002-08-15, 14:28:06 


23.9 


2.35x10^1 


2.50x10^1 



'^relative to the HEG range: 1.5-15 A 
^relative to the MEG range: 2-24 A 
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The HETG spectra and the corresponding X-ray hghtcurves for the observations can be 
found in TDP04; we also discuss there the variabihty observed for some of the stars. The 
source sample covers very different stellar characteristics, providing a wide view of X-ray 
emitting stellar coronae; e.g., the X-ray luminosities (in the hetgs bandpasses; see TDP04) 
span more than five orders of magnitude, from the relatively weak emission of Proxima 
Centauri with a few 10^^ erg s~^, up to the very high luminosity (~ 6 x 10^^ erg s~^) of the 
giant HD 223460. 



4. Analysis 

The data used here were obtained from the Chandra Data Archiv^ and have been 
reprocessed using standard CIAO v3.2.1 tools and analysis threads. Effective areas were cal- 
culated using standard CIAO procedures, which include an appropriate observation-specific 
correction for the time-dependent ACTS contamination layer. Positive and negative spectral 
orders were summed, keeping heg and meg spectra separate. For sources observed in sev- 
eral different segments, we combined the different observations and analysed the coadded 
spectra. For IM Peg, as noted above, a number of observations are available with which to 
probe optical depth properties at different times and orbital phases; IM Peg will be discussed 
in more detail in §5.11 

Spectra were analyzed with the PINTofALsi IDL§ software jKashvap fc Drakel l200o[ ) 



using the technique of spectral fitting described in TDP04. We measured the spectral line 
intensities of the Nex Lja {2p ^-P3/2,i/2 — Is and Ly/3 {3p ^^3/2,1/2 — l-s ^'S'1/2) transitions 

from both HEG and MEG spectra, and the O VIII lines from MEG spectra alone (since the 
latter lie outside the heg wavelength range). 

In order to take into account the mild dependence of the Lya/Ly/3 ratio on plasma 
temperature, we also measured the intensities of the Neix and OVH resonance line, r 
{ls2p ^Pi — Is^ ^5*0 ), providing us with an estimate of a representative temperature through 
the Lya/r ratio. There are two potential problems with this approach. Firstly, where reso- 
nant scattering is relevant, both Lja and r transitions can be depleted, and in such a case 
the Lya/r ratio might deviate from its expected theoretical behaviour. Secondly, the coro- 
nae in which these X-ray lines are formed are expected to be characterised by continuous 



''http:/ /cxc. harvard.edu/cda 
'■'http://hea-www.harvard.edu/PINTofALE 
^Interactive Data Language, Research Systems Inc. 
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ranges of temperatures, rather than by a single, isothermal plasma. However, both of these 
concerns prove unproblematic. 

Figure [2] illustrates the theoretical temperature dependence of the Lya/r ratio and of 
the Ly/J/Lya ratio for Ne lines (O lines show analogous behaviour). The Lya/r is much more 
temperature-sensitive than Ly/?/Lya, and the relatively small deviations from the optically- 
thin case that might be expected in stellar coronae can only incur small errors in temperature 
that will have a negligible impact on the predicted Ly/3/Lya; ratio. 
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Fig. 2. — Temperature sensitivity (in the isothermal case) of the Lya/r ratio compared to the 
Ly/3/Lya ratio, for neon hnes. The corresponding oxygen hnes show analogous behaviour. 
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The accurate temperature determination from the Lya/r ratio strictly holds only for 
isothermal plasma, whereas stellar coronae are characterised by a thermal distribution of the 
plasma, so that this diagnostic would give us a temperature weighted by the emission mea- 
sure distribution (DEM). By computing line ratios for non- isothermal plasma models, and 
interpreting as if isothermal, we can estimate the errors incurred by an isothermal assump- 
tion. We considered two different sets of DEMs: those derived from actual observations 
of some of the stars in our sample, and simple DEM models in which the e mission mea 



sure is proportional to T^/^ (as expected for simple hydrostatic loop models; iRosner et al. 



Scelsi et al. 



1978 ). or proportional to T^^^ (as observed in some stars such as 31 Coni, 
2004 . and reproduced by some hydrodynamic loop models, iTesta et al.ll2005l ). For these 
models we used peak temperatures varying from 10^'^ K to 10^'^ K. Observed DEMs were 



culled from the literature for the following sources: AB Dor (ISanz-Forcada et al.l l2003bl ) 



HD 223460 JXesta et al.lboOTh. 31 Com. 3 G et, Vel ( 



Z CrB, e UMa (ISanz-Forcada et al.ll2003ar). 44 Boo flBrickhouse fc Dupree 



Garcia-Alvarez et al 



20061 



199 



ER Vul, 
, UX Ari 



(ISanz-Forcada et al.ll2002i), II Peg (IHuenemoerder et al.l 120011). A And 



Sanz-Forcada et al, 



2OO2I ). AR Lac (IHuenemoerder et al 



20031 ). HR 1099 (iDrake et al.ll200ll ) 



Given the line ratio computed for each DEM model, we inverted the isothermal relation 
to obtain T(Lya/r). Using that T, we then obtained the theoretical isothermal Lya/Ly/? 
ratio (Lya/Ly/3 [T(Lya/r)]) and compared to the synthetic ratio (Lya/Ly/5 [DEM]). We 
find that the isothermal assumption is a good predictor of the ratio for both theoretical 
DEMs and for representative stellar models: Figure [31 shows that the isothermal and DEM 
ratios are in very good agreement, with differences between the two generally amounting to 
a few percent. We conclude that the Lya/r temperatures are quite adequate to assess the 
appropriate expected Lya/Ly/? ratio, even if the real plasma temperature distributions of 
the coronae in our study are far from isothermality. 
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Fig. 3.— Le/t ^)ane/.- DEM models (DEM oc T^/^ T^/^) ^gg^ ^gg^ ^j^g validity of the 
isothermal approximation. Right panel: The Nex Lya;/Ly/5 ratios obtained for model (left 
panel) and observed DEMs (see text; black empty symbols) compared with the Lya/Ly/? 
ratios expected for isothermal temperatures diagnosed from the DEM Ne Lya/r ratios. 
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4.1. Deblending of Nex Lya and O VIII Ly/5 

Ne X Lya and O Vlll Ly/3 lines are affected by blending of iron lines unresolved at the 
HETGS resolution level. Specifically, significant blending is expected for O Vlll Ly/? by an 
Fe XVIII line {2s^2p^{^P)3s ^Pg/s -25^2^5 ^Py2, A = 16.004A), and, to a lesser extent, for the 
Nex Lja line by a nearby Fexvil transition {2s'^2p^{^P)4:d ^Pi - 2s^2f ^5*0, A = 12.124A). 



O VIII Ly/3 — In order to estimate the often significant contribution of Fe xvill to the 
O VIII Ly/5 spectral featur e, we used the same rn ethod described in Paper I, which was also 
subsequently adopted by iNess fc SchmittI (120051 ) in an analysis of spectra of the classical 
T Tauri star TW Hya. We estimated the intensity of the 16.004A Fexvill line by scaling 
the observed intensity of the slightly stronger neighbouring unblended Fexvill 16.071 A 



(2s22p4(3p)3s ^P5/2 - 25^2/ 2pg^2) transition. 1Gu| fl2003b[ ) has recently pointed out the 



similar behaviour of Fe L-shell lines originating from 3s and 3p upper levels in respect 
to the indirect excitation processes of radiative recombination, dielectronic recombination, 
and resonance excitation. The 16.004 A and 16.071 A transitions originate from similar 3s 
upper levels and their ratio should therefo re not deviate gr eatly from current theoretical 
predictions. The APED (v. 1.3.1) database (jSmith et al.ll200ll ) lists their theoretical ratio as 
0.76 at the temperature of the Fexvill population peak (~ 6.8 MK), with a decrease of only 
a few percent up to 10 MK or so; this range covers the expected formation temperatures of 
Fe XVIII in our target stars. 

We investigated this scaling factor empirically by examining the departure of the re- 
sulting deblended Lya/Ly/? ratios from the theoretical value as a function of the ratio of the 
observed Fexvill 16.071A to Ly/5 line strengths. If the deblending scaling factor is correct, 
there should be no residual slope in the resulting data points; a positive slope would instead 
indicate an over-correction for the blend and a negative slope an under-correction. We found 
that a slope of zero was obtained for a scaling factor of ~ 0.70 — in good agreement with the 
APED value within expected uncertainties. This is illustrated in Figure |H 
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Fig. 4. — Ovill Lya/Ly/? ratios (relative to the theoretical ratio) before (left) and after 
(right) the correction for the contamination of the Fexvill line blending with the O Ly/? 
hne, plotted vs. the relative intensity of the Fexvill hne at ~16.07 A and the O Ly/3 be- 
fore the deblending process. The effect of the blending Fexvill line is clear from the left 
plot. The right plot shows the effectiveness of the deblending procedure which ehminates 
the strong correlation between the departure from the theoretical Lyct/Ly/? value and the 
Fexvill/OLy/? ratio before the deblending. 
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Ne X Lya — In analogy with the procedure used for the O Vlll Ly/3 we searched for isolated 
and strong Fexvil lines to estimate the amount of blending from Fexvil A = 12.124A in 
the Nex Lya line. A potentially good candidate would be the nearby 12.266A Fexvil line 
(2s^2|}^(^P)4(i ^Di — 2s^2p^ ^•S'o), which has a similar intensity and behavior as a function 
of temperature; however this transition is rather weak in most of the observed spectra, and 
can be affected itself by blending with a close Fexxi transition (12.284A) barely resolved by 

HETG. 

We therefore chose to use a larger set of Fexvil lines including stronger transitions. 
The intensity of the Fexvil line blending with the Nex Lya is estimated by scaling the 
observed intensity of ot her Fexvil lines (12 .266, 15. 014 15.261 16.780, 17.051, 1 7.096A) by 
the ratio expected from lLandi fc Gd (120061 ) (see also lGull2003al : iLandi et al.ll2006l ). assuming 
logT = 6.8. These calculations include indirect processes involving the neighboring charge 
states, and are quite successful at reproducing the relative intensity of these Fexvil lines. 
Previous calculations (including those adopted in APED v. 1.3.1) fail here, with, e.g., the 
ratio of the 15A line to the 17A lines being overestimated by about a factor 2. The detailed 
comparison of the observed Fexvil line ratios with the predictions of different theoretical 
calculations are addressed in a paper in preparation. 

The measured line fluxes after application of deblending corrections, together with the 
statistical errors, are listed in Table [2l 



5. Results 

The observed Lja/Ljj3 ratios and the temperature derived from the Lya/r ratios are 
listed in Table [3] and illustrated in Figure O Measurements obtained from the meg spectra 
generally result in smaller statistical errors due to the higher S/N; however, we find a general 
agreement of heg and meg result^. The main differences of the present work with respect to 
the analysis of Paper I, where results were presented for 4 sources (II Peg, IM Peg, HR 1099, 
and AR Lac), are that the line fiuxes have been remeasured in the reprocessed data, Nex 



^The issue of some systematic discrepancies between HEG and MEG measurements has been addressed 
briefly in TDP04, and it is thoroughly discussed in http://space.mit.edu/ASC/calib/heg_meg/, where sys- 
tematic discrepancies of 8-10% are found in the 10-12A range (meg fluxes being lower than corresponding 
HEG fluxes at the same wavelength) therefore not significantly affecting the Lya/Ly/3 ratios. We note that 
our measured HEG fluxes are in agreement with meg measurements within 2a for the large majority of cases; 
the only significant exceptions are the Nex Lya of TZ CrB and of HR 1099 (3, 5a respectively) whose HEG 
fluxes are ~ 10% larger than the corresponding meg fluxes in agreement with the comparisons presented in 
http://space.mit.edu/ASC/calib/heg_nieg/. 
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Lya fluxes have been corrected for blending, and the deblending procedure for O Vlll Ly/3 has 
been shghtly refined (see previous section). Also, IM Peg fluxes listed here refer to the total 
(~ 192 ks) HETGS spectrum including ObsID 2530-2534 not yet publicly available when we 
carried out the analysis presented in Paper I (where we used the first three 25 ks pointings; 
see §5. II for a detailed discussion). 



-17- 



Table 2. Measured line fluxes (in 10 ^ photons cm ^ s 



Source 






Ne 








O 




Fe 






Lya 


a 




Ly^J 


r 


Lya 


Ly/3 


r 


Fe xvii"^ 


Fe XVIII 




12.132A 


10.239A 


13.447A 


18.967A 


16.006A 


2I.6O2A 


12.124A 


I6.O71A 




HEG 


MBG 


HEG 


MEG 


MEG 


MEG 


MEG 


MEG 


MEG 


MEG 


AU Mic 


374 ± 50 


356 ± 15 


51 ± 10 


44 ± 4 


208 ± 15 


990 ± 40 


1,34 ± 20 


243 ± 40 


21 ± 5 


67 ± 11 


Prox Cen 




69 ± 20 




6.5 ± 2.9 


47 ± 9 


306 ± 28 


38 ± 14 


100 ± 70 


6 ± 2 


19 ± 13 


BV Lac 


278 ± 22 


262 ± 10 


30 ± 6 


32 ± 3 


209 ± 8 


875 ± 19 


137 ± 11 


330 ± 30 


24 ± 5 


73 ± 9 


AB Dor 


668 ± 50 


682 ± 24 


96 ± 16 


86 ± 9 


332 ± 27 


1400 ± 50 


200 ± 22 


280 ± 50 


41 ± 9 


168 ± 16 


TW Hya 


97 ±19 


95 ± 8 


12 ± 8 


9.6 ± 2.8 


155 ± 14 


266 ± 30 


38 ± 12 


114 ± 40 


6.9 ± 2 


11 ± 8 


HD 223460 


168 ± 20 


137 ± 11 


22 ± 6 


16.0 ± 2.7 


20 ± 7 


229 ± 22 


27± 11 


48 ± 24 


10 ± 2 


42 ± 10 


HD 111812 




58 ± 5 




7.9 ± 1.6 


14 ± 3 


124 ± 16 


16 ± 7 


40 ± 30 


8.8 ± 2 


48 ± 5 


/3 Get 


419 ± 24 


408 ± 10 


53 ± 10 


48 ± 4 


146 ± 10 


764 ± 29 


62 ± 23 


140 ± 40 


116 ± 21 


426 ± 21 


HD 45348 




54 ± 4 




4.7 ± 2.3 


35 ± 7 


127 ± 22 


16 ± 8 


40 ± 30 


11.6 ± 3 


30 ± 5 


^l Vel 


173 ± 21 


159 ± 15 


16 ± 13 


18 ± 7 


95 ± 7 


400 ± 30 


45 ± 16 


46 ± 27 


62 ± 12 


163 ± 10 


Algol 


898 ± 50 


832 ± 40 


136 ± 23 


120 ± 8 


196 ± 20 


1318 ± 50 


166 ± 30 


190 ± 60 


63 ± 14 


268 ± 23 


BR Vul 


218 ± 22 


201 ± 12 


26 ± 8 


24.3 ± 2.8 


74 ± 7 


362 ± 29 


39 ± 12 


66 ± 40 


24 ± 5 


88 ± 9 


44 Boo 


572 ± 60 


637 ± 24 


76 ± 11 


68 ± 8 


254 ± 20 


1450 ± 60 


122 ± 22 


266 ± 60 


53 ± 11 


161 ± 16 


TZ CrB 


1126 ± 30 


1043 ± 24 


146 ± 14 


127 ± 6 


516 ± 20 


2110 ± 50 


178 ± 26 


300 ± 50 


131 ± 25 


445 ± 20 


UX Ari 


730 ± 40 


767 ± 20 


126 ± 17 


107 ± 8 


245 ± 18 


926 ± 50 


147 ± 21 


158 ± 40 


14.5 ± 3 


55 ± 13 


5 XJMa 


419 ± 27 


394 ± 13 


42 ± 8 


47 ± 4 


309 ± 16 


1510 ± 60 


l(i9 ± 22 


49(i ± GO 


83 ± 17 


185 ± 15 


II Peg 


1300 ± 70 


1200 ± 30 


183 ± 24 


191 ± 10 


356 ± 20 


1950 ± 70 


350 ± 30 


250 ± 50 


18 ± 4 


86 ± 16 


A And 


667 ± 30 


666 ± 14 


80 ± 11 


84 ± 6 


189 ± 14 


984 ± 40 


113 ± 17 


110 ± 30 


24.6 ± 6 


137 ± 12 


TY Pyx 


341 ± 29 


312 ± 6 


46 ± 12 


40 ± 5 


142 ± 15 


468 ± 50 


51 ± 21 


96 ± 50 


25 ± 5 


106 ± 17 


AR Lac 


648 ± 40 


621 ± 17 


104 ± 15 


94 ± 6 


182 ± 20 


888 ± 50 


94 ± 23 


130 ± 50 


34 ± 7 


129 ± 14 


HR 1099 


1960 ± 40 


1730 ± 18 


255 ± 14 


230 ± 7 


566 ± 18 


2680 ± 50 


367 ± 21 


410 ± 40 


58 ± 12 


246 ± 14 


IM Peg 


408 ± 25 


383 ± 10 


62 ± 7 


60 ± 3 


79 ± 8 


490 ± 30 


75 ± 12 


50 ± 20 


10.6 ± 2.5 


45 ± 8 



^Ne LycK fluxes after the deblending with the Fe XVII line at 12.124A. 
''O Ly/5 fluxes after the deblending with the Fe xvill line at 16.004A. 

'^Fexvii line (at 12.124A) flux estimated by scaling the measured fluxes of the other Fexvii lines. 
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Fig. 5. — Lya/Ly/?, both from HEG (filled symbols) and from MEG (empty symbols) vs. 
temperature, as derived from the ratio of the Lya and the He-like resonance line. The 
solid line mark the theoretical ratio from APED, expected for isothermal plasma at the 
corresponding temperature. As indicated on the right plot, different symbols are used for 
different classes of sources. Lya/Ly/^ ratios from HEG are shifted in T by 2% to separate 
them from the MEG measurements. 
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Tablc 3. Lya/Ly.5 photon ratios and plasma temperature, with la errors. 



Source Ne O 





rjp a 


Lya/Ly/3 


A [a 


•] 

J 


y a 


Lya/Ly/3 


AH t 

L J 




[IQS K] 


HEG 


MEG 


HEG 


MEG 


[10^ K] 


MEG 


MEG 


AU Mic 




7.3 ± 1.6 


8.1 ±0.8 


-0.4 


0.1 




7.4 ± 1.1 


-0.7 


Prox Ccn 






10.6 ± 6 




0.4 


. . _ _|_n 7 
'^•^-0.9 


8.1 ± 3.1 


-0.1 


EV Lac 


. ri-l-O 1 


9.3 ± 1.9 


8.2 ± 0.8 


0.5 


-0.1 


0.0_Q 2 


6.4 ±0.5 


-4 


AB Dor 


5 

^" ' —0.2 


7.0 ± 1.3 


7.9 ±0.9 


-0.6 


0.3 


^■^-0.4 


7.0 ± 0.8 


-0.9 


TW Hya 


o q+0.1 

''•^-O.l 


8.1 ±5.6 


9.9 ±3.0 


-0.1 


0.3 


o n+0.3 
o.U_o.4 


7.1 ±2.6 


-0.7 


HD 223460 


q 1+1-3 


7.2 ±2.1 


8.6 ± 1.6 


0.3 


1.2 




8.5 ±3.6 


0.2 


HD 111812 


7 9+0.7 
' -^-0.6 




7.3 ± 1.6 




0.1 


3.5^0;^ 


8.3 ±4.0 





p Get 


6.ot°:? 


7.9 ± 1.5 


8.5 ±0.7 


0.2 


1.2 




12.3 ±4.2 


1.0 


HD 45348 






11.2 ±5.4 




0.6 


o c+o.g 

'J-'J-l.S 


7.9 ±2.3 


-0.2 


/i Vel 


4 S+0.2 


10.8 ±6.8 


8.8 ±3.5 


0.4 


0.1 




8.9 ±3.2 


0.5 


Algol 




6.6 ± 1.2 


6.9 ±0.6 


-0.4 


-0.3 


A 8+0-3 


7.9 ± 1.5 


0.2 


ER Vul 




8.4 ±2.7 


8.3 ± 1.1 


0.3 


0.8 


4.4+0-8 


9.3 ±3.0 


0.5 


44 Boo 


^- ' -0.2 


7.6 ± 1.4 


7.9 ± 1.0 


-0.1 


0.2 




11.9±2.2 


1.9 


TZ CrB 


t; c+0.2 
"3-3-0.1 


7.7 ±0.8 


8.2 ±0.4 


-0.1 


0.9 




11.9 ± 1.8 


2.2 


UX Ari 


6.6t°;? 


5.8 ±0.9 


7.1 ±0.6 


-1.7 


-0.3 


4.4t°1 


6.3 ± 1.0 


-1.3 


^ UMa 




10.0 ±2.0 


8.4 ±0.8 


0.8 


0.0 


o t-+0.1 
0.0_Q 2 


8.9 ± 1.2 


0.4 


II Peg 


6.9t°;? 


7.1 ±1.0 


6.3 ±0.4 


-0.2 


-2.8 


A fi+0.4 
*-"-0.4 


5.6 ±0.4 


-5 


A And 


"J-"-0.3 


7.0 ± 1.0 


6.6 ±0.4 


-0.4 


-1.8 




8.7 ± 1.4 


0.9 


TY Pyx 


r 7+0.3 

"3- ' -0.2 


7.4 ± 1.8 


7.8 ±0.9 


-0.1 


0.1 


4 2+0-6 


9.2 ±3.9 


0.4 


AR Lac 




6.2 ±1.1 


6.6 ±0.6 


-1.0 


-1.2 


4 6+0-6 

+ •"-0.8 


9.4 ±2.4 


0.8 


HR 1099 




7.7 ±0.4 


7.5 ±0.2 


0.7 


0.5 


4.6±°J 


7.3 ±0.4 


-0.6 


IM Peg 




6.6 ±0.8 


6.4 ±0.3 


-0.5 


-1.8 


5.5t?;« 


6.5 ± 1.1 


-0.8 



*T derived from the Lya/r line ratio diagnostics. 

^Discrepancy between Lya/Ly/3 measured and theoretical value in units of a. 
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The observed Nex Lya/Ly/5 ratios follow closely the APED theoretical predictions as 
a function of temperature. The values for the O Vlll ratios (right panel of Fig. [5]) show 
slightly more scatter, with some departures below theory and some l-2a excursions above. 
We note that an enhancement of the Lya/Ly/5 ratio is possible for the particular geometry 
in which an emitting region is optically-thick in some directions but is optically-thin in the 
line-of-sight of the observer (see also the radiative transfer study of iKerr et al.ll2004l ). Such 
a situation is disfavoured by the isotropic nature of the scattering process in which any line 
enhancements through scattering are diluted, roughly speaking, by the ratio of solid angles 
of optically-thin to optically-thick lines-of-sight. In the case of loop geometries, such a ratio 
is much larger than unity. We therefore interpret these small upward 1-2 a excursions as 
normal statistical fluctuations. 

In the case of the Nex lines, the meg Lja/LjjS ratio for II Peg is ~ 3cr lower than 
the expected value; the ratio derived from heg is consistent with the meg measurement 
but does not depart significantly from theory. The only sources showing OviII Lya/Ly/5 
significantly (> 3a) lower than the theoretical ratio are II Peg and EV Lac. We note that 
the low Lya/Ly/? ratios for IM Peg discussed in Paper I were found in the spectra obtained 
in the first three pointings while here we analyze the whole set of observations; in §5. II below 
we discuss the variability observed in Lya/Ly/? ratios for IM Peg. 

One other interesting case is the T Tauri star, T W Hya, that shows very high density. 



Tie ^ 10 cm~ , at the temperature of O Vll lines (e.g.. lKastner et al.ll2002l ). In this case, the 
high densities are generally believed to arise from an accretion shock, rather than coronal 
loops; significant scattering effects could in principle provide interesting constraints on the 
dimensions of the shock region. Unfortunately, for this source the relatively low S/N results 
in large error bars in the Lya/Ly/5 ratios th at do not provide any useful constraints (as was 
also noted in the analysis of Fexvil lines by lNess fc Schmittll2005l ). 



5.1. IM Peg 



The different Chandra hetg observations of the RS CVn system IM Peg provide an 
opportunity for studying the optical depth and its possible v ariation with time . This system 
has been analyzed in great detail at optical wavelengths by iBerdyugina et al.l (Il999l . 120001 ). 
These works depict a scenario with stellar spots concentrated mai nly close to the polar 



regions, similar to those found for many other active systems (e.g., ISchuessler et al.l Il996 



Hatzes et al.lll996l : IVogt et al.lll999l : iHussain et al.ll2002l : iBerdvurina et al.l 119981 ). 
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Chandra observed IM Peg eight times over ~ 2 orbital periods (see Tab. [T]). The X-ray 
hghtcurve obtained from this sequence is illustrated in Figure [6] and shows significant modu- 
lation that is possibly related to the or bital and rotational peri od (these being essentially the 
same; Prot = 24.39 d, Porb = 24.65 d, IStrassmeier et al.lll993l ). In Paper I we analyzed the 
first three 25 ks hetgs observations of IM Peg, publicly available at that time, and we found 
remarkably small photon path lengths based on optical depth effects. The compact emission 
regions implied might be associated with active regions revealed by optical Doppler imaging 
studies. Any significant net line photon loss to resonant scattering out of the line-of-sight 
would be expected to be sensitive to the orientations of the emitting structures. Modulation 
of the Lya/Ly/? ratio with time could provide important new structure and morphology 
diagnostics. Such regions might also be expected to change on relatively short timescales as 
a result of flaring activity and consequent re-alignment of their defining magnetic fields. 
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Table 4. Lya/Ly/S photon ratios with 1 a errors. 



Obs. ID 


Nex 


VIII 


HEG 


MEG 


MEG 


2527+2528 


5.3 + 0.8 


5.5 + 0.5 


5.1 + 1.2 


2529+2530+2531 


6.9+1.7 


6.5 + 0.6 


7.5 + 3.0 


2532+2533+2534 


6.2+1.2 


7.1 + 0.6 


7.0 + 2.1 



10 



32 



0.0 



X 



10 



31 



0.5 



Phase 
1.0 



1.5 



2.0 




100 200 300 400 
Time (lO^sec) from MJD 52456.669 



Fig. 6. — Lightcurve of IM Peg Chandra observations; both HEG {filled squares) and MEG 
{empty squares) Lx are presented. HEG values are shifted by 10^ s on the time axis with 
respect to the corresponding meg measurements. The phases corresponding to the absolute 
times are indicated at the top of the plot. 
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4 F 




10.0 10.5 11.0 11.5 12.0 

wavelength [A] 




10.0 10.5 11.0 11.5 12.0 

wavelength [A] 



Fig. 7. — Nex Lya,/3 spectral region for the three chosen time segments of the observation 
of the RS CVn system IM Peg, as indicated in the plots. Both heg [upper panel) and 
MEG [lower panel) spectra are presented, normalized to the intensity of the LjjS line at 
10.24 A. For better readability the spectra corresponding to the second (2529+2530+2531) 
and the third (2532+2533+2534) segment are shifted in wavelength with respect to the first 
(2527+2528) by +0.03 A and +0.06 A respectively. Note the pronounced difference in Lya 
and Ly/5 relative intensity in the three spectra especially between the first and the third 
spectra. 
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Here we present an analysis of the complete dataset, exploring possible temporal vari- 
ability of the optical depth of the coronal plasma. Unfortunately, the single IM Peg exposures 
provided counts sufficient for reliable intensity measurement only for the strongest lines in 
the spectrum, such as Nex Lya. In order to probe secular change in the Lya/Ly/3 ratio, 
we therefore measured lines from spectra coadded from contiguous sets of two or three ob- 
servations that sample different phases of the stellar period: ObsIDs 2527 and 2528; 2529, 
2530, and 2531; and 2532, 2533, and 2534. The spectra in the Ne and O Lja,P regions for 
these three portions of the observation are shown in Figure [7] and [Hi The measured Lja/Lj(3 
ratios are hsted in Table HI 
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16 17 18 19 

wovelength [A] 



Fig. 8. — Ovill Lya,/3 spectral region of meg spectra of IM Peg for the three chosen 
time segment of the observation, in the same format of the plots in Figure [7l The shift in 
wavelength of second and third spectra is of +0.06 A with respect to the preceding spectrum. 
A clear trend in the Lya/Ly/5 line ratio is present. 
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The Lya;/Ly/5 ratios, while statistically all consistent with one another, are suggestive 
of different conditions in the different portions of the observation. Both Ne and O Lya /Ly/5 
ratios in heg and meg spectra are lower than the corresponding APED theoretical value 
in the first segment, as already discussed in Paper I and noted earlier. The rest of the 
observation is characterized by Lya/Ly/3 values compatible with theory. The measured 
ratios are illustrated as a function of the hardness ratio {HR = {H — S)/{H + S), where 
H is the flux integrated in the 2-9A wavelength band, and S is the flux integrated in the 
9-25A band) in the different segments of the observation in Figure O There is no obvious 
correlation of Lya/Ly/? with spectral hardness that might suggest, e.g., that the former 
changes as a result of plasma temperature changes. These results show that it might be 
worthwhile to explore, when the data quality allows it, temporal variability of optical depth 
properties in stellar coronae, expected to some extent on the basis of their dependence on 
the line-of-sight and coronal geometry other than on the physical conditions of the emitting 
plasma. 
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■ 1 ■ 1 ■ ■ 1 ■ 


■ 1 ■ 1 


- Al 

- ^) 

NeX 

-- OVIII 

1 .... 1 . . 


A N 

A N 
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1 .... 1 . 


eX (HEG) - 
eX (MEG) 
IVIII (MEG) " 



+ 

X 



LD 



II 

X 



0.08 
0.06 

0.04 

0.02 

0.00 

■0.02 

■0.04 
■0.06 



H 
S 



2.0A-9.0A 
9.0A-25.0A 



[] 



0.0 



0.5 



1.5 



2.0 



1.0 
Phase 

Fig. 9. — LyQ;/Ly/3 ratios measured from IM Peg spectra for the three time segments selected 
for the analysis (2528+2528, 2529+2530+2531, and 2532+2533+2534; upper pane!), together 
■with the hardness ratio, HR {lower panel). In the upper panel the horizontal lines represent 
the APED values expected for the Lja/Ljf3 ratio of O {dashed line) for a temperature of 
5 MK, and of Ne {dotted line) for a temperature of 9 MK. For better readability Lya/Ly/? 
ratios of Nex from HEG, and of O vill are shifted ■with respect to meg Nex Lya/Ly/^ ratios 
by +0.05 and +0.1 respectively on the phase axis. 
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5.2. Path length estimates 



The data presenting more reliable evidence for optical thickness are the spectra of IM Peg 
from the first two observations, where Lya of both oxygen and neon appear to be depleted 
with respect to the corresponding Ly/3 transition (deviation of ~ 1.8a and ~ 3a respec- 
tively), and the spectrum of II Peg characterized by an Ovill Lja/Ljj3 ratio ~ 5a lower 
than the theoretical value. As pointed out in Paper I, the discrepancies between observed 
and theoretical Lya/Ly/5 line strength ratios in IM Peg and II Peg cannot be explained 
by photo electric absorption along the line-of-si ght on the basis of their measured H column 
densities (IMewe et al.lll997l : iMitrou et al.lll997l ). It is worth discussing whether these transi- 
tions, and in particular their ratios, might be affected by non-equilibrium conditions that are 
possibly relevant for these active stars undergoing frequent fiaring activity. Non-equilibrium 
effects might be responsible for cha nges in the O yil (j ratio ((/ + ^)/t) observed for EV Lac 
between quiescent and flare phases (jMitra-Kraevll2006l ). It is well-known that the G ratio is 
sensitive to deviations from coronal equilibrium, because re combination pr ocesses contribute 



signiflcantly to / and i but are almost negligible for r (e.g. lPradhanlll985l ). Even though re- 
combination processes have non-negligible effects on Lya and Ly/? (10-15%), they contribute 
to a very similar extent to the two transitions, and therefore their ratio is not expected to 
change signiflcantly under mildly non-equilibrium conditions. Furthermore, we note that for 
the typical conditions of the coronal plasma in these active stars the expected times cales 
of these effects are only of the order of hundreds of seconds (e.g., iGolub et al.lll989l ). i.e. 
very short with respect to our integration times of several tens of kiloseconds. We therefore 
interpreted the discrepant ratios in terms of a relative depletion of the Lya line flux due to 
resonance scattering processes. 

Another source apparently showing signiflcant effects of resonant scattering is EV Lac, 
whose O VIII Lya/Ly/3 ratio is more than 3a lower than t he corresponding t heoretical value. 
This result is supported to some extent by the flndings of iNess et al.l (120031 ): in their survey 
of coronal optical depth properties based on the analysis of Fe xvil transitions, EV Lac is 
the only source deviating signiflcantly from the typical values found for all other coronae. 
However, iNess et al.l (120031 ) did not consider this result robust due to discrepancies between 
HEG and MEG measurements. 

For the stars with Lja/Ljj3 departing from the theoretical va lues we can derive an 
estimate of the photon path len gth using the iKaastra fc Mewe ( 1995 ) approximation to the 
escape probability formalism of Kastner fc Kastnei ( 199ol ). as described in § 3.1 of Paper I. 
The observed Lya/Ly/? can be expressed in terms of the line center optical depths r'^, f- values 
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Z*^, and photon pathlength as 



{I Lya / 1 Ly 13) ohs _ 1 
{iLya/ lLyl3)th 3 



+ 



l + 0.43C(£)/f 1 + 0.43C(£)/f /2 



1 + 0.43C(£)/^ 



(2) 



where C{i) ■ = and the hne center optical depth is given by 



1 ifi in-i4 "-i /I "-H w 

1.16 ■ 10 Az — \f\ —He 



(3) 



where n-Jnei is the ion fraction (from iMazzotta et aLlll998l ). Az is the element abundance, 
n-R/ue ~ 0.85, / is the oscillator strength, M the atomic weight, T the temperature, and He 



the electron density. The Lya and Ly/3 f-values are 
Morton|[2003l ). 



0.2776, and /f = 0.05274 (e.g. 



We use the electron densities, Ue, derived from the diagnostics of the He-like triplets 
TDP04. The coronal abundances assumed for II P eg and IM Peg are d iscussed in Paper I; 
for EV Lac we assume O/H=8.40, as derived by iFavata et al.l (l2000bl ). expressed on the 
usual spectroscopic logarithmic scale in which X/H= log[n(X)/n(H)] + 12, where n(X) is the 
number density of element X. 

In Table [5] we list the values obtained for the path length estimates, ir, and for com- 
parison, we list the stellar rad ii and loop lengths expected for a standard hydrostatic loop 
model (e.g. 
and densities. 



Rosner et al.lll978l . RTV hereafter) corresponding to the observed temperatures 
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Table 5. Path length derived from measured Lya/Ly/? ratios. 



Source 








-^RTV ^ 






Es " 








[cm] 


[cm] 






[erg cm~^ s~^] 


IM Peg 


Ovm 




1.5 • IQi" 


2.2 • 10^ 


0.017 


- 0.006 


- 7.2- 10^ 




Nex 


[HEG] 


2.1 • 10« 


2.8 • 10^ 


0.00024 


- 0.0003 


- 5.2- 10^° 




Nex 


[MEG] 


1.7-10*^ 


2.8 • 10^ 


0.00019 


- 0.0004 


- 5.3-101° 


II Peg 


VIII 




9.5 • 10^ 


1-109 


0.04 


- 0.021 


- 1.2- 10* 


EV Lac 


VIII 




1.6- 10^ 


2.6- 10* 


0.058 


~ 0.016 


~ 1.9- 10* 



''Loop length from RTV scaling laws: Lrtv ~ T3/[(1.4 x lO^)^ ■ p]. 
'^Path length as fraction of the stellar radius. 

^Coronal surface filling factor, defined as f, = A/A^ = {V/tr)IA^ = [EM / {n%)]l A^. 
•^Estimate of surface heating flux (see text for details). 
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6. Discussion 

This work present a detailed and extensive study of tlie O Vlll and Nex Lya/Ly/? ratios 
in a sample of active stars, paying careful attention to the presence of blends and to the 
validity of theoretical line ratio predictions. Our analysis shows that optical depth effects 
are generally negligible in the disk-integrated X-ray spectra emitted by stellar c oronae over 



a wide range of activity, i n line with previous studies based on Fexvil lines (e.g.. iNess et al. 



20031 : lAudard et al.l l2004j ) . We argued in ^ that the coronal abundance patterns exhibited 
by most active stars renders the Ne and O Lyman lines more sensitive diagnostics of optical 
depth than the Fe lines, and in this regard our study extends the results of previous surveys. 

However, our study has also yielded convincing indication of line-of-sight photon loss 
through resonant scattering in the spectra of three stars, IM Peg, II Peg and EV Lac. The 
detection of significant X-ray optical depth is particularly interesting, because it provides us 
with insights into the characteristic dimensions of the emitting coronal structures. 

Before proceeding, we note here that a simple escape probability analysis does not 
include the scattering source term itself, and so photons scattered into the line of sight that 
could enhance the observed line strength are not included. For instance, in the case of a 
spherically symmetric corona in which the line optical depth were significant, scattering into 
and out of the line of s ight would be balanced and hne strengths not affected (see e.g.. 



Wood fc Raymond! I2OOOI ). Strictly, then, the photon path, £, entering into the equations in 



5. 21 should be interpreted as a lo wer limit to the true photon path length. This is potentially 



important. iPhillips et al.l (l200ll ). for example, used the lack of appreciable depletion of the 
Fexvil resonant line at ~ IS.OlA to deduce an upper limit of 3000 km for size of emitting 
regions in the corona of Capella; however, such a measurement in reality only provides an 
upper limit to the lower limit, therefore not providing a constraint. Indeed, the lack of 
evidence for line quenching through scattering in most of our sample stars does not imply 
that scattering is not a significant source term, but merely that the scattering geometry is 
such that there is no net loss or gain of photons in the line-of-sight. 

Instead, a positive detection of resonance scattering has interesting implications and 
points to a non-uniform "aspect ratio" of the dominant coronal emitting regions: for any net 
line depletion to occur, an emitting structure must generally be more elongated along the line 
of sight than in the perpendicular direction. In the context of a corona comprising plasma 
contained by magnetic loops, there are different ways in which this can be interpreted: (1) 
scattering loss results from the structure and viewing angle of the loops themselves; (2) 
scattering loss arises because of a particular conglomeration of loops viewed at a particular 
angle. Each case leads to fundamentally different interpretations of observed scattering. In 
the former case, a loop with random orientation must generally be viewed from the top, in 
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which direction the photon path length through the plasma is largest. Photon loss through 
resonance scattering then implies that loops are preferentially placed on the stellar surface 
facing the observer (a perfect alignment of loops seen edge-on on the stellar limb could 
produce a similar effect, though such a chance alignment is unlikely). In the latter case, the 
coronal structures or active regions must be placed preferentially on the stellar limb. 



6.1. Scattering within individual loops 

With the proviso that a scattering-derived path length is a lower limit to the true 
emitting region size, we note that the photon path lengths implied by the analysis of line 
ratios ( §5.21) observed in II Peg, IM Peg and EV Lac and listed in Table [5] are all much 
smaller than the corresponding stellar radii, but are also about an order of magnitude larger 
than expected from RTV loops. 

Under the assumption that scattering within individual loops arises when they are 
viewed from the top, we can interpret the path length estimates, ir, in terms of the coronal 
scale height. With knowledge of the total emitting volume we can also derive an estimate of 
the surface filling factors. Emitting volumes can be estimated using the emission measures 
implied by the different lines, V = EM/nl, combined with electron densities, Ue, derived 
in TDP04. The surface filling factors, fs, are then given by fg = Aj A^, = {V/ir)/A^^. The 
derived filling factors, listed in Table [5], are very small, and especially so for the hotter plasma. 
They are also an order of magnitude smaller than those previously derived in TDP04 based 
on RTV loops — a direct consequence of the optical depth scale height being commensurately 
larger than those suggested by simple quasi-static loop models. 

The small filling factors we find here have interesting implications for the surface heating 
flux requirements. A very rough estimate of the heating required to sustain the observed 
physical conditions of the plasma can be obtained from the RTV scaling laws. While our 
estimated loop lengths are significantly longer than suggested by RTV models, these relations 
should still suffice for the purposes of estimation. The volumetric heating per unit time, E, 
IS given by E ~ 10^ • p^/^L-^/^, where p is the plasma pressure and L the loop length. By 
assuming L = £r we find for the the surface flux values of order of 10^ erg cm^^ s^^ for 
O VIII lines, and several 10^° erg cm~^ s~^ implied by the Nex lines in IM Peg (Table |5]). 
These compare with typic al surface heating rates f or the cores of solar active regions of a 



few 10^ erg cm ^ s ^ (e.g.. IWithbroe &: Noved 119771 ). 



In the comparison of ir with Lrtv we assumed to some extent that ir is a reasonable 
estimate for the coronal scale height (i.e. for the loop length). However, an alternative 
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interpretation is possible in which the actual coronal loops are larger and ^r represents only 
an estimate of the length of the part of the loop containing plasma at the characteristic 
temperature of the quenched lines. In such a scenario, the values found for 1^ can suggest 
the presence of cooler loops with maximum temperature around 3 MK, similar to solar loops, 
coexistent with loops with much higher maximum temperature (>10 MK). In these hotter 
loops, the plasma emitting the Nex lines (T > 6-7 MK) would be confined in the lower 
portion of the loop, which is characterized by higher density. In order to estimate the region 
of the loop occupied by the plasma emitting the observed line. A/, we can use the equations 
for the conductive flux: 

and from RTV relations for a radiative power loss function P{T): 

FJAI ~ nlP{T). (5) 

Assuming T is the temperature of maximum formation of the line k, T^.^^, and AT is the 
width in temperature of the line emissivity curve (~ 0.3 dex), we can solve the equations for 
A/: 

A/ ~ 5 ■ lO^cm for O Vlll 

A/ ~ 2 ■ 10 W for Nex 

The resulting values of A/ are rather close to L/jtv (see Table [5]), obtained assuming T^^^ 
as maximum temperature of the loop, and are still much smaller than ir- We can conclude 
then, that the values derived for the path length do not seem to agree with the hypothesis 
of standard uniformly heated quasi-static loop models. 



6.2. Scattering within active regions 

The maximum path length, /max, through a spherically-symmetric corona of height h on 
a star with radius i?^ and surface filling factor fg (assuming this filling factor does not vary 
significantly with height) is 

/max = 2fJh^ + 2hR,. (6) 

Considering typical parameters found for stellar coronae, we can estimate from the above 
equation whether or not we expect significant photon scattering in stellar coronae, regardless 
of whether we see a net photon loss. From Eq. [3] we can estimate the optical depth at the limb 
for Ne and O Lya lines. Assuming the temperature of maximum formation of the line, i.e. 
about 6 MK and 3 MK for Ne and O respectively, and typical density of about 10^'^ cm~^ and 
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10" respectively (see e.g., TDP04. lNess et al.ll2004h . we derive r(Ne) ~ 5.2 x lO-^^-Zmax 
and r(0) ~ 3.6 x 10-^° ■ /^ax- 

The height derived for coronal structures with different techniques (see ^Hfor references) 
is < on the Sun a typical coronal height is closer to i?^/10. Assuming h = and 

R-^ = Rq we obtain /max ~ fs ■ ~ /s • 6 x 10^° cm. 

Under these assumptions, the values of the optical depth at the limb, for a spherically 
symmetric corona, are therefore r(Ne) ~ 30/s and t(0) ~ 20/^ (or r(Ne) ~ 80/^ and 
r(0) ~ 55/s if we assume h = Ri,/2, and R^ = Rq). For the typical filling factors lower 
than a few percent (e.g., TDP04) we obtain optical depth r < 1 both for Ne and for O. 



6.3. Correlation of scattering with stellar and coronal parameters 

While only three of our stellar sample present a significant case for resonance scattering, 
it is possible that trends of the departure from the theoretical Lya/Ly/3 ratio with funda- 
mental stellar parameters can be found. We have examined the departures of observed from 
theoretical Lya /Lj/S ratios as a function of X-ray surface flux, surface flux in the Lya line, 
Lx/Lhoi, filling factors, and plasma density. These comparisons are suggestive of correla- 
tions between quenching of Lya photons and both Lx/i^boi and the density- sensitive ratio 
of strengths of the forbidden and intercombination lines, f /i, of Mgxi. 

These correlations are illustrated in Figure [TO] where we show the measured to theoretical 
ratios of Ovill Lya/Ly/? of our sample as a function ofLx/i^boi (top) and Mgxif/i {bottom). 
Error-weighted linear fits to the data in these figures yield slopes of —0.11 ± 0.04 and 0.25 ± 
0.06, respectively. The sources presenting evidence of optical depth effects are the stars 
characterised by the highest activity level (Lx/i^boi) and the highest plasma densities (i.e. 
lowest f/i) in our sample. Such correlations are what are expected based on Eqn. [31 Optical 
depth is proportional to the product of electron density and typical path length within an 
emitting region, n^i. For a given fixed volume emission measure, nlV ~ n^l^, the optical 
depth varies as nl^^ or and so increases with increasing plasma density. In the case 

of Lx/I^boi, an increase in Lx can arise through either i or ng, such that any increase in Lx 
might typically be expected to lead to greater scattering optical depth. While only a small 
handful of our measurements present truly significant detections of resonance scattering, the 
evidence for trends of increasing r with Lx/L\^o\ and Mgxi //i as is expected adds confidence 
to the interpretation of the Lya/Ly/? ratios in these terms. 
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Fig. 10. — Ratios of measured to theoretical O Vlll Lya/Ly/? ratios plotted vs. Lx/I>boi {top) 
and vs. the MgXI forbidden to intercombination line ratio (from TDP04; bottom), which is 
a diagnostic of plasma density (low f/i ratios correspond to high density). 
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7. Conclusions 

We have investigated the optical thickness of stellar coronae through the analysis of 
Lya and Ljj3 lines of hydrogen-like oxygen and neon ions, in Chandm-UETG spectra of a 
large sample of active stars. Our study indicates that most stellar coronae are characterised 
by negligible visible signs of optical depth, in agreement with the results of previous studies 
based on Fexvil hnes. This indicates that coronae are either in general optically-thin, or 
that for cases in which optical depths reach of order unity or higher the geometry does not 
strongly favour lines-of-sight showing net Lya photon loss. 

We do find evidence of significant optical depth in the O Vlll Lyman lines of the RS CVn 
binary II Peg, and of the single M dwarf EV Lac; the RS CVn binary IM Peg also shows 
depletion of both Ne and O Lya/Ly/? ratios as compared with theoretical predictions, and 

our analysis indicates that it is a transient effect present only in part of the observations. 
The detection of significant optical depth allows to derive an estimate for the photon path 
length and therefore for the typical height of the corona. The size of coronal structures 
derived for all three sources is of the order of a few percent of the stellar radius at most, 
implying very small coronal filling factors and high surface heating fiuxes. We searched for 
correlation with basic stellar parameters and coronal properties and we find that the sources 
presenting evidence of significant optical depth are at the high end of activity level, with 
Lx/i^boi at the saturation limit, and high densities in their hot plasma, as revealed by the 
Mg XI He-hke triplet hnes. 

For the stellar sample as a whole, we also find evidence of increasing quenching of 
Lya relative to Ly/3 as function of both Lx/ L^oi and the density-sensitive Mg xi forbidden 
to intercombination line ratio. Such a trend is expected in the scenario in which optical 
depths are significant but generally small: viewing geometry rarely favours large net photon 
enhacements, but for favourable lines-of-sight photon depletion is expected to increase with 
both increasing Lx and increasing plasma density. 
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